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Abstract 

Using data collected with the CLEO detector operating at the CESR e'^e" collider at i/s = 3.97- 

4.26 GcV; wc investigate 15 charmonium decay modes of the ^'(4040), ip{4160), and y(4260) reso- 
nances. Wc confirm, at 11a significance, the BaBar y(4260) tt~^ tt~ J / ip discovery, make the first 
observation of y(4260) -> n^ir^J/ip {5.1a), and find the first evidence for y(4260) K+R-J/ip 
{3.7a). We measure e+e~ cross-sections at -^i = 4.26 GeV as (j(7r+7r~ J/'i/')=58l}Q±4pb, 
(T(7r°7r°J/^)=23l!l8^iblpb, and a{K'^K~J/il^)=9'^^±lph, in which the uncertainties are statis- 
tical and systematic, respectively. Upper limits are placed on other decay rates from all three 
resonances. 
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The region at center-of-mass energies above charmonium open-flavor production threshold 
is of great interest to theory due to its richness of cc states, the properties of which are not 
well-understood. Prominent structures in the hadronic cross-section are the ip{3770), the 
?/'(4040), and the ^^(4160) Q. Their main characteristics are large total widths, two orders 
of magnitude larger than for the lower- lying cc states of J^'" = 1 , and weaker couplings 
to leptons than the J/ip and ip{2S). Decays to closed-charm final states are not favored due 
to the availability of open-charm channels. 

Recently, observations of new charmonium-like states in the same energy region have been 
reported j^. Most of these have been observed to decay through open charm. However, an 
enhancement in the invariant mass spectrum of the closed-charm tt^tt^J/iIj final state has 
also been observed by BaBar in initial state radiation (ISR), e"''e~ — > 'y(7!''^Tc~J/ip) j^. A 
weaker signal was observed in the decay B K{7t^7t^ J/ip) \^. The observed lineshape can 
be described by a single resonance, termed y(4260), of mass M = 4259 ± 8tl MeV, width 
Ftot = 88 ± 23t^MeV, and coupling Tee x i3(F(4260) ^ n+n-J/ip) = 5.5 ± 1.0lg:?eV. It 
is located quite unexpectedly at a local minimum of the hadronic cross-section. Since it 
is observed in ISR, the new state must have J"^*" = 1 , and therefore it can be studied 
directly in e"'"e~ collisions at threshold. No other evidence for a resonance at this mass has 
been identified j^, leaving the existence and possible charmonium-like nature of this state 
uncertain. Many interpretations have been suggested; to be compatible with the absence 
of a corresponding enhancement in open charm production, most favor an unconventional 
explanation of y(4260), such as hybrid charmonium td,0|; tetraquarks or hadronic 

molecules jlllllj, C)ne proposal argues that F(4260) is conventional: It identifies 
y(4260) with the 4'{4:S) vector cc state, and relies upon interference effects to produce the 
dip in open-charm cross-section and a hypothesized large coupling to tt~^tt^J/iP of the ipi'^S), 
commonly associated with the '?/'(4040). 

To further clarify the nature of y(4260), investigation of both open and closed charm is 
necessary. Here, we report production cross-sections measurements of 16 final states contain- 
ing a J/ip, ip{2S), Xcj, or 0, in the ?/'(4040), -0(4160), and F(4260) energy region, motivated 
by the range of experimental tests suggested so far 

II 0, i i H Q Ha, Q El- we 

use data taken during a scan of center-of-mass energies ^/s =3.97-4.26 GeV, complement- 
ing a sample of an integrated luminosity / Cdt = 281 pb~^ at v/s = 3.773 GeV previously 
acquired. Collisions were registered with the CLEO detector [l5| at the CESR e"*"e~ col- 
lider The scan data naturally separate into three regions: the ?/'(4040) {^/s = 3.97- 
4.06 GeV, fCdt = 20.7pb-^), the V^(4160) (4.12-4.20 GeV, 26.3 pb"^), and y/s = 4.26 GeV 
(13.2 pb~^). Fig. d'a) shows a profile of / Cdt vs. ^/s and Fig. [Tfb) the Born-level Breit- 
Wigner lineshapes l| for the four resonances between y^=3.7 and 4.4 GeV, also indicating 
the grouping of scan points. We estimate the number of resonances produced by folding 
together the luminosities with the resonance Breit-Wigner cross-sections jl^, including ra- 
diative corrections [M, lij, and arrive at (93±ll)xl0=^ [(115±15)xl03] '0(4040) [0(4160)] 
mesons, where the dominant errors are the uncertainties on resonance parameters. 

The CLEO detector features a solid angle coverage of 93% for charged and neutral 
particles. The charged particle tracking system operates in a 1.0 T magnetic field parallel 
to the beam axis and achieves a momentum resolution of ~0.6% at momenta of IGeV/c. 
The Csl crystal calorimeter attains photon energy resolutions of 2.2% for E^=l GeV and 5% 
at 100 MeV. Particle identification is performed with the specific ionization loss [dE/dx) 
and the Ring Imaging Cherenkov detector (RICH). Muons of momentum p > 1 GeV are 
separated from pions by their penetration of the calorimeter, solenoid coil, and up to three 
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36-cm-thick slabs of magnet iron for subsequent detection by wire chambers behind each 
slab. The integrated luminosity was measured using e"'"e~, 77, and events l2ll|. 

The final states analyzed here are listed in Table HI We require all particles in each final 
state to be reconstructed, and four-momentum conservation is enforced. Mass windows for 
the following light hadron decays are set, based on MC studies: 7r° 77 (110-150 Me V); 
V ^ 11, V ^ 7r+7r-7r° (450-650 Me V); r]' 7r+7r-?7, 7]' 'yn+n- (930-980 Me V); u 
Tr+vr-vrO (730-830 MeV); and ^ K+R- (1.00-1.04 GeV). 

We identify a. J/ip or ip{2S) through its decay into £+£~=e+e^ or ■ A lepton 

candidate is identified by the ratio of the energy deposited in the calorimeter, to the 
measured momentum, p. Muon pair candidates must satisfy E/p <0.25 for at least one of 
the tracks and E/p <0.5 for the other; both electrons must have E/p >0.85 and a dE/dx 
consistent with the value expected for an electron. Pions faking muons are additionally 
suppressed by requiring at least one muon candidate per muon pair to leave a signature 
in the muon system. A lepton pair is classified as a J/ip [ip{2S)] — > candidate if the 
invariant mass of the decay products lies within 3.04-3.14 [3.64-3.73] GeV. We also use the 
decay ipi^S) vr+vr" J/?/', where we require Miji^n' J /ip) =3.64-3.73 GeV. 

A 7Xci,2 11 J /'4' decay is tagged by the highest energy photon in the event having an 
energy within (-60,-^40) MeV of E^={s - M\^^)/2^s. Events with M(77) in the 7r° or 77 
mass ranges are excluded. Any calorimeter shower other than the two radiative transition 
photon candidates must have E <50MeV. We select 7r"'"7r~7r°Xcj, Xcj — ^ iJ/'^ events 
(J = 1,2) by demanding that M{'yJ/'il}) as well as the mass recoiling against the 7r^7r~7r° 
match My.^j |1] within ±20 MeV. 

For K^K~ J /ip, the kaons have momenta that are too low for use of the RICH detector, so 
TT^ rejection is achieved by requiring that at least one of the candidates have momentum 
in the range 0.2-0.5 GeV/c and have dE/dx within three standard deviations of the expected 
value for a K^. For 7r+7r~0 and Xco 2{tt^it^) or tt^tt^K^K^ , candidates must 
be positively identified as a kaon using a likelihood based upon both dE/dx and RICH 
responses, but none of the vr^ candidates in either decay can be so identified as kaons. For 
the oJXco niode, both the mass of the XfEs decay products and the mass recoiling against the 
oj must lie within 50 MeV of M(xco) 111- Similarly, for ir^n'cj), the mass recoiling against 
7r+7r~ must lie in the range 0.94-1.10 GeV. 

For -K^Ti^J/il) and 7r+7r^^/'(2S'), we require that M(7r+7r^) > 350 MeV and that neither 
pion candidate be identified as an electron via E/p and dE/dx as above. These cuts suppress 
e^e~ — >■ £"'"£"7, 7 — >■ e^e~ events in which the e^e~ pair from the photon conversion is 
mistaken for the tt+tt". For n^J/ip as well as the rfJ/ip and iXci,2 modes that end in 
77e"'"e~, background from Bhabha events is diminished by requiring cos^e+ < 0.3. 

For modes with a J/ip or 'ip{2S) we restrict the missing momentum k, computed from the 
measured event momenta according to Eq. (6) of Ref. [23|, to further reduce background. 
Signal events will have /c ~ 0; selection windows are set separately for each mode accord- 
ing to the measurement resolution predicted by MC simulation, and range from ±10 MeV 
{K+K^J/i)) and ±15 MeV (tt+tt-J/V^) to (-35,±75) MeV {-K^-K^J/i)). 

Backgrounds in the XJ/ip modes from e^e~ 7*— s>hadrons are estimated from the 
yields seen in the data with the restrictions on M{[^[~) and k adjusted to correspond 
to masses either below or above Mj/^. The sideband windows, \k\ < 150 MeV and ei- 
ther M(£"'"£~)=2.90±0.15 GeV or 3.30±0.15 GeV, are considerably wider than the signal 
regions, allowing for the accumulation of more statistics in the background estimates, and 
therefore are scaled down by the sideband/signal window width ratios prior to subtrac- 
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tion from the signal yields. The sideband scale factors are verified in MC simulations of 
e"'"e~ — >■ 7* — >-hadrons. We perform a mass sideband subtraction for the Xco, as a result 
of which all observed events can be attributed to background, and for and r]', also re- 
ducing the yield substantially; photon energy sidebands are used to estimate background 
in 7Xcj- Background in rnrJ/tp, riJ/ip, and 7Xcj modes produced in '^ip{2S) decay (with 
M{'il)(2S)) ~ s/s ^ 3686 MeV) is summed along with that from the sidebands in Table HI 
and comprises 10-40% of the total background in those modes. Other modes could be sub- 
jected to similar subtractions, but due to the absence of significant event populations, we 
forgo this option, and use uncorrected yields to compute (conservative) upper limits. 



The EVTGEN event generator [2^, which includes final state radiation [23[, and a geant- 
based ji^ detector simulation are used to model the physics processes. The generator 
implements a relative S-wave (P-wave) configuration between the vrvr [j] or 7r°) and the 
J/ if) or %1){2S). Detection efficiencies (Tabled) range from 4-38%, not including the effects 
of the intermediate branching fractions 

TT+TT" J/'?/', Xcj iJ 1'^-, or Xco — >hadrons, but including those for 7r°, r], uj, r]' , and decay. 
Already included are the effects of ISR, which reduce efficiencies by relative fractions of 
8-21%. 

The radiative return process e^e~ — 'jipi'^S) —^XJ/ip results in final states which are 
nearly identical to some of our signal modes, and thereby affords an opportunity to verify 
our understanding of efficiencies, background, and luminosity. To gather such events, we 
alter only the /c-windows of the event selections, as such events will congregate not near 
= but rather around = {s — M^^2S))/'^V^- The cross-section for this process can be 
calculated by integrating the convolution of a Breit-Wigner lineshape (approximated by a 
5-function) with the ISR kernel W from Eq. (28) of Ref. |3: 

^(e+e- ^ 7V^(25)) = x Wis, 2k,/ ). (1) 

This curve as well as the cross-section measurements using il){2S) — > XJjjh, X=7r+7r^, 7r°7r°, 
and 7/ jl^, are shown in Figure [Tfc). The CLEO value for Vee[ip{2S)] j22] sets the scale of 
the theoretical curve in the figure. The predicted number of observed radiative return events 
for the three channels together in all the scan data is 820±21, which compares favorably to 
825±29 events seen. 

Table E] lists cross-section results for the /c = region. Cross-section central values and 
errors are shown when the statistical significance (the likelihood that the observed event 
yield is due entirely to background) exceeds 2. Sex; otherwise upper limits at 90% confidence 
level (CL) are shown. Poisson fiuctuations of the background are taken into account in the 
computation of statistical uncertainties. Cross-sections for e^e~ irnJ/ip vs. y/s are shown 
in Fig. ^d); selected missing momentum and dipion mass distributions appear in Figs. |21 
and El respectively. The yields for n^n^cj) are consistent with the rate observed [1^ from 
e+e" ^ 7* at ^^=3.67 GeV. 

Systematic uncertainties arise from the following sources: luminosity (2%), charged par- 
ticle tracking (1% per track), particle identification (1% per high momentum e, tt, or K; 5% 
per K in K^K^ J/tp), 3% per fi pair for muon chamber modeling, and a mode-dependent 
contribution amounting to 50% of the total estimated background. Statistical uncertainties 
dominate. The total systematic error excluding the background uncertainty amounts to ap- 
proximately 10% for most channels (up to 30% for Xcj modes and 7T~^7i~(j)). The only charmo- 
nium channels with more than 2.5a statistical significance are, at -\/i=4.260 GeV, 7r"^7r~ J/V' 
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(11a), Ti^Ti^J/i> (5.1a), and K+R-J/ip (3.7a); in the ^^(4160) dataset, n+n-J/ip (3.6a) and 
n^^^J/tp (2.6a); and in the V^(4040) dataset, vr+vr-J/V' (3.3a). The ^(4160)-region yields 
of TT^n~J/ip and n^ir^J/ilj are consistent with being due entirely to the y(4260) low-side 
tail j^. No compelling evidence is found for any other decays in the three resonance regions, 
and corresponding upper hmits on cross-sections (and, for ip^AOAO) and ■?/'(4160) datasets, 
on branching fractions) are set. In particular, we find i3('?/'(4040) -k^-k^ J /ip) <0.4% and 
i3('?/^(4160) — i> 7r^7T~J/ip) <0.4%. These correspond to partial widths of less than 0.4 MeV in 
both cases, to be compared with the central values for ip{2S) and ?/'(3770) of 0.10 MeV and 
0.44 MeV |l|, respectively. While statistics are low, no prominent narrow features emerge in 
M{7r~^TT~), and the distribution is somewhat softer than the '0(25')-like MC prediction. 
This analysis provides a high-significance confirmation of the BaBar signal of 7r^Tc~J/tp. 



The observation of the n^ir^J/ip mode disfavors the XcjP^ molecular model ll|. The fact 
that the Ti^ir'^J/ip rate is about half that of 7!-~^n~J/ip disagrees with the prediction of the 
baryonium model Our evidence of significant K~ J /ip production is not compatible 
with these two models either. No evidence of a large h'^tt' J/ijj signal from the '?/'(4040) 
is observed, making the conventional Y (4:260)=^ (AS) assignment [lj| less attractive. The 
results are compatible with hybrid charmonium interpretations 0,01; but open-charm studies 
will be required to make more definitive conclusions. A large coupling to Di(2430)°Do and 
a small one to D^Dg signals hybrid charmonium 9, whereas a dominant D^Dg (DD) could 
favor a esse 0,9 {cqqc (sj) tetraquark model. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the A. P. Sloan Foundation, 
the National Science Foundation, and the U.S. Department of Energy. 
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TABLE I: For each mode e^e~ X, for three center-of-mass regions: the detection efficiency, e; 
the number of signal [background] events in data, Ng [A^b]; the cross-section a{e'^e~ — X); and 
the branching fraction of ■0(4040) or '0(4160) to X, B. Upper Umits are at 90% CL. '-' indicates 
that the channel is kinematically or experimentally inaccessible. 
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FIG. 1: (a) / Cdt vs. ^/s (circles), (b) Born-level Breit-Wigner cross-sections for resonances in 
this energy region; the y(4260) curve has an arbitrary vertical scale, (c) e~^e~ — > ^ip{2S) cross- 
section measurements vs. ^/s from 'k'^'k~ J /ijj (circles), J /il) (squares, dashed lines), and 77J/V' 
(triangles) overlaid with the theoretical prediction as described in the text, (d) Cross-sections for 
e'^e~ — > ■n'^TT~ J/ip (circles) and J/ijj (squares, dashed lines) vs. ^/s. Some points in (c) and 
(d) are offset in ^/s by lOMeV for display purposes. 
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FIG. 2: The missing momentum [k] distribution for tt~^tt~J/iJj (top), Tr^TT^J/ip (middle), and 
K~^K~ J/ip (bottom) in the data at y/s = 4.26 GeV (circles), and the signal shape as predicted by 
MC simulation (solid line histogram) scaled to the net signal size. 
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FIG. 3: The efficiency-corrected dipion invariant mass distribution for the vr+vr" J/t/' final state in 
the signal region at ^/s = 4.26 GeV for the data (circles) and the signal shape (solid line histogram) 
as predicted by '0(2S')-like MC simulation scaled to the net signal size. 
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